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SUMMARY 

This report illustrates the use of NASTRAN® modal synthesis capability 
for some small examples. A classical truss problem is examined and the 
results for accuracy are compared to existing results from other metnods . 
fins proDiem is examinea using botn fixed interface modes and tree interface 
modes. Tne solution is carried out tor an applied dynamic load down as far 
as recovery of forces in individual members as a function of time. Another 
small beam problem is used to compare different means of "combining" 
substructures . 

INTRODUCTION 

During the past twenty years, a body of technology has developed within 
the general field of structural dynamics that has been identified by the 
term modal synthesis. Modal synthesis is a Rayleigh-Ritz approach using 
systematically derived displacement functions. It is used to formulate and 
solve the large eigen problems which arise in dynamic analysis of complex 
structural systems. Solutions are approximate in the sense that the motion 
of the structure is constrained to linear combinations of a limited nunber 
of modes or displacement functions characterizing the behavior of indepen- 


dent substructures. 


Several researchers have formulated various modal synthesis procedures 
in an attempt to reduce computation errors and minimize computer costs. 

Hurty developed the first modal synthesis method capable of analyzing 
structures with redundant interface connections in references 1 and 2. He 
treated the structure as an assembly of connected components, or substruc- 
tures, each of which is analyzed separately to derive a set of modes or dis- 
placement shapes from which a set of generalized coordinates applicable to 
the complete structure is synthesized. Craig and Bampton (ref. 3) simpli- 
fied Hurty' s formulation by combining two groups of coordinate functions 
which Hurty had defined separately. A number of survey papers have been 
written by Hou, Goldman, Benfield and Hruda in references 4 to 7. Some 
methods are found to be more suitable for certain applications than others. 
Yet, experience has shown that no single approach is generally preferred 
over the others. 

The complexity of aerospace structures increased enormously during the 
last two decades. A new challenge is presented by the proposed space sta- 
tion (ref. 8) in that it is an evolving structure that cannot be ground 
tested because final configuration may not be known when the first component 
is put into space. Therefore, the component mode synthesis method may be 
applied for the dynamic analysis of such large structure system in space. A 
widely used tool for structural analysis, the NASTRAN® computer program, 
contains a modal synthesis capability but, other than the demonstration 
problem presented in reference 5, little is publicly known about its 
capabilities. 

The purpose of the present report is to examine some of the capabili- 
ties of this program. This is done by examining two simple problems, a 
truss and a beam. 

NUMERICAL EXAMPLES 

The modal synthesis procedure in NASTRAN® is applied to two simple 
structures. One is a redundant truss confined to lie in a plane but free to 
move in this plane. It is composed entirely of ROD elements (no bending 
stiffness for all) . This example is used to examine convergence character- 
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istics of the modal synthesis procedure and also to illustrate the transient 
response capability all the way down to obtaining stresses in rod members 
as a function of time. The second example is a free-free beam. It is used 
to examine different ways to "combine" substructures to yield frequency for 
the total structure. 

Truss Example 

The redundant truss example is the one used in reference 5 to compare 
eight different modal synthesis procedures. The full truss model is shown 
in figure 1(a) and its two components shown in figure 1(b). Component A 
consists of five equal bays and has a total of 18 joints. Component B con- 
sists of four equal bays and has a total of 15 joints. All members in the 
components have identical properties. At the interface of the components in 
the full truss model, the vertical member has twice the area of other mem- 
bers. Basic geometric and material properties are presented in table I 
along with the prescribed load for a transient response analysis. An addi- 
tional run was made with the full model subdivided into three components 
with three bays in each component. 

The basic run sequence and substructure operation are shown in figure 
2. In the figure capitalized letters inside of rectangular blocks indicate 
names of psuedostructures used in the analysis. Capitalized letters adja- 
cent to, or on, the flow diagram indicate the names of modules that perform 
a certain function in the computer program. At the top of figure 2, the 
Phase 1 operations formulate the finite element stiffness and mass matrices 
using Rigid Format 2. For the convergence study the Phase 2 runs on Rigid 
Format 3 were repeated using a different number of modes from the individual 
components. Also Phase 2 runs were using free interface modes as well as 
the fixed interface modes. A limited amount of data is presented for three 
components and naturally a Phase 1 run must be made for this component. 

A transient response analysis was made on this free-free truss struc- 
ture for an axial load applied to the right end of the truss. The load was 
applied for 0.12 seconds and then removed. In order to apply a load at grid 
point 42 in component B, this grid point must be included on a BOUNDARY 



card. Thus, additional degrees of freedom are created corresponding to this 
point. The structure was represented by eight modes from component A, six 
modes from component B, and the eight interface modes for a total of twenty- 
two modes. The modes for the individual component were determined with the 
interface fixed. The standard procedure will obtain displacements back in 
the individual component. However, member forces and stresses are not de- 
termined automatically, but can be obtained through a simple procedure in a 
few steps. In the first step a run is made with DIAG 17 turned on to put 
the EMAP sequence on the punch file with an EXIT scheduled after statement 
1. A small substructure deck is included to allow the appropriate commands 
that interface to the Substructures Operating File (SOF) to be generated. 
This punch file is subsequently saved and altered to replace the RECOVER 
module with the SDR2 module which can recover element forces and stresses. 
The listing of this DMAP sequence and run stream is contained in Appen- 
dix A. 


Beam Example 

This example consists of a beam composed of seven components as shown 
in figure 3(a) . All subbeams have a constant length, area and uniform mass 
properties. Each component consists of ten equal elements and has a total 
of 11 joints as shown in figure 3(b) . Basic geometric and material 
properties for each subbeam are presented in table II. A lumped mass 
formulation is used (no rotary inertia) and, therefore, there are 213 
stiffness degrees of freedom in the problem, but only 142 eigenvalues. 

Three different ways of "combining" substructures are illustrated in 
figures 3(c), 3(d), and 3(e). The basic run sequences and substructure 
operations for each case are shown in figures 4 thru 6. For all cases, the 
sub structuring Phase 1 operations formulate the finite element stiffness and 
mass matrices for subbeam A using Rigid Format 3. The structural matrices 
contained in BBASIC, CBASIC, ..., FBASIC are generated as needed by using 
EQUIV operation. The basic subbeams are reduced to modal coordinates and 
combined together following the procedures shown in figures 4 thru 6. The 
eigenvalues of the total beam are obtained by using the MRECOVER command. 

The driver decks and sample bulk data for cases 1, 2 and 3 are listed in 
Appendices B, C and D. Only fixed interface modes were used but two sets of 
runs were made using a different number of modes from the subbeams. 
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RESULTS 


For assessing the accuracy of the modal synthesis procedure, two and 
three truss components with fixed or free interface connection are run to 
determine frequencies and compared to results for full model. Percentage 
errors in frequency for the combined systems of 12, 20, 28 and 36 degrees of 
freedom are shown in tables III thru VI. Here degrees of freedom include 
not only the number of flexible modes used but also any interface modes. 
Thus, for example, for 12 degrees of freedom results, since there are six 
interface modes, only six flexible modes can be shown. Based on the lowest 
frequency criterion then four modes were chosen from component A and two 
modes from component B. 

Figures 7 thru 11 are nondimens ional plots that indicate the relative 
accuracy obtained by modal synthesis procedures. Also shown on the figures 
are results taken directly from reference 5 in which several other proce- 
dures are compared. From figures 7 to 10 it can be seen that modes derived 
with the interfaces fixed yield better results than modes derived with the 
interface free. 

For the transient response run the percentage error in displacement for 
grid points 41, 42, and 43 of component B are shown in table VII. These 
results were produced from the 20 degrees of freedom model. The axial force 
in elements 111-113 and 143 of component B are shown in table VIII. 

The full beam shown in figure 2 was run to determine its natural fre- 
quencies and used as a comparison of results obtained with the various "com- 
bination" procedures. Table IX shows the percentage error in frequency for 
the various "combination" procedures when 62 degrees of freedom are used. 
These 62 degrees of freedom correspond to approximately 47% of the total 
degrees of freedom in the full model. All three "combination" procedures 
yield good results. However, case 1 uses considerably less CYBER 75 CPU 
time than the other two cases (53.8 CPU seconds corresponds to 65.3 seconds, 
59.1 seconds, respectively). Another run for case 1 was made using 19% of 
total degrees of freedom, and 55% frequencies were obtained with less than 
1% error in frequency. 
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Typical frame width 
(see fig. 1(b)) 


a = 1.015 m (40") 


Typical frame height 
( see fig. 1(b)) 


h = 0.762 m (30") 


Cross-sectional area of 
members 


A = 1.935 cm 2 (0.3 in 2 ) 


Young's modulus 


E = 1.422xl0 4 K| (10 7 psi) 


Density 


P = 272.517 Kg-sec 


(2 . 5x1 0 -4 lbf-sec 2 ) 


Transient loads 


P 42 = 2. 2x1 0 3 Kg(10 3 lbf) 0< t<0 . 1 2S 


0 


t>0 . 1 2S 











Table II. Beam Geometric and Material Properties 

I 
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Typical component length 
[ see fig. 3(b)) 

1 = 2.54 m (100") 

Cross section of beam 

A = 3.613 cm 2 (0.56 in 2 ) 

Young's modulus 

E = 1.422xl0 4 Kg (10 7 psi) 

7 

Density 

P = 282.437 Kg-sec^ 

m 

(2.591xl0 -4 lbf-sec 2 
in 4 

Total beam length 

L = 15. 78 m (700") 








for TVo Modal Synthesis Models Using 12 Degrees of Freedom 



Mode No. 


Full Truss 
(Hz) 

Free Interface 
(%) 

Fixed Interface 
(%) 

65.7771 

29.88 

0.0015 

136.3306 

4.32 

0.0022 

175.5505 

21.9 

0.0314 

202.7780 

7.35 

0.0198 

260.3387 

3.49 

0.0536 

316.2614 

5.12 

0.0227 

334.1522 

61.41 

4.21 

347.1668 

142.05 

6.439 

388.1286 

ii 

183.78 

0.97 






































Table IV. Frequency for Full Truss and Percent Error in Frequency for Ttoo 
Modal Synthesis Models Using 20 Degrees of Freedom 


Mode 

No. 

Ful 1 Truss 
(Hz) 

Free Interface 
(%) 

Fixed Interface 
(%) 

3 Components 
Fixed Interface 

1 

65.7771 

19.23 

0.00074 

0.00351 

2 

136.3306 

2.82 

0.00044 

-0.03425 

3 

175.5505 

8.2 

0.0087 

-0.00915 

4 

202.7780 

2.67 

0.0087 

0.02232 

5 

260.3387 

2.23 

0.0091 

-0.00355 

6 

316.2614 

2.0 

0.0078 

-0.03588 

7 

334.1522 

0.65 

0.75 

0.00521 

8 

347.1668 

3.9 

0.088 

0.00469 

9 

388.1286 

0.3 

0.23 

-0.01105 

10 

394.1834 

0.3 

0.1 

-0 . 00029 

11 

414.9853 

1.9 

0.18 

-0.00924 

12 

451.2226 

8.57 

0.078 

-0.00182 

13 

466.3475 

8.5 

0.14 

0.00130 

14 

504.7402 

7.8 

0.41 

0.01524 

15 

507.2363 

39.7 

1.32 

0.03394 

16 

537.3632 

58.4 

2.0 

0.01038 

17 

575.3048 

114.65 

0.7 

0.00005 
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Table VI. Frequency for Full Truss and Percent Error in Frequency for Two Modal Synthesis Models 
Using 36 Degrees of Freedom 
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Table VII. Transient Response and Percent Error in Displacement 









































































































Times 


111 



0.012 


137.4126 













Table IX. Percent Frequency Error Using 62 Degrees of Freedom 
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Table IX. Percent Frequency Error Using 62 Degrees of Freedom (concluded) 
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Figure 1. Truss Model 
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Runs 1 and 2 
Phase 1, RF2 



Run 3 

Phase 2, RF3 


Run 4 

Phase 3, RF9 



Run 5 
DMAP 



Figure 2. Substructure Formulation Tree and Solution Sequence 
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(a) Total beam model. 
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2.54 m 
( 100 ") 


(1) | (2) | (3) | (4) | (5) 1(6) | (7) | (8) | (9) | (10) 
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(b) Representative finite element model of any component. 
i 2.54m | 

p-(100")-**| 

i a i a i S I ~0" I E I ■ F — 'I E 


(c) Case 1 - All components combined simultaneously. 
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(d) Case 2 - Components combined sequentially. 
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(e) Case 3 - Components combined in pairs. Pairs then combined sequentially. 


Figure 3. Total Beam Model and Various Subdivided Representations. 
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Figure 5. Case 2 Subbeam Formulation Tree and Solution Sequenc 
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Figure 6. Case 3 Subbeam Formulation Tree and Solution Sequence 




















Figure 7. Comparison of Methods with Frequency Error of 








Figure 8. Comparison of Methods with Frequency Error of 0.5% 
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Figure 9. Comparison of Methods with Frequency Error of 
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Figure 10. Comparison of Methods with Frequency Error of 5.0%. 
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Figure 11 . Comparison of Methods with Frequency Error of 10 . 0 %. 









APPENDICES 


APPENDIX A. Driver decks and sample bulk data for two components 
truss problem. 


nastran files = umf $ coc and tbm 

ID DEM2031 uNASTRAN 

APR DISP.SU8S 
SOL 2.0 
TIME 3 

CEND 

SUBSTRUCTURE PHASE 1 
PASSWORD = MDLSYN 

S0F(1) = FT19.5Q0* NEW S CDC AND IBM 
NAME = A0ASIC 
SOFPRINT TOC 
ENOSUBS 

TITLE = TRUSS DYNAMIC ANALYSIS USING AUTOMATED MOOAL 
LA8|L= SUBSTRUCTURE 1, RUN 1, PHASE 1, RF 2 
SUBTITLE = NASTRAN DEMONSTRATION PROBLEM NO. 2-3-1 
BEGIN BULK 

2 
3 

12 
13 
22 
23 

32 

33 

42 

43 

52 

53 
11 
12 
13 
21 
22 
23 

31 

32 

33 

41 

42 

43 

51 

52 

53 
13 
21 
11 
23 
31 
33 
41 
43 
51 
53 


CROD 

1 

1 

1 

CROO 

2 

1 

2 

CROD 

11 

1 

11 

CROD 

12 

1 

12 

CROD 

21 

1 

21 

CROD 

22 

1 

22 

CROO 

31 

1 

31 

CROD 

32 

1 

32 

CROD 

41 

1 

41 

CROD 

42 

1 

42 

CROD 

51 

1 

51 

CROO 

52 

1 

52 

CROD 

111 

1 

1 

CROD 

112 

1 

2 

CROD 

113 

1 

3 

CROD 

121 

1 

11 

CROD 

122 

1 

12 

CROD 

123 

1 

13 

CROO 

131 

1 

21 

CROD 

132 

1 

22 

CROD 

133 

1 

23 

CROD 

141 

1 

31 

CROD 

142 

1 

32 

CROD 

143 

1 

33 

CROD 

151 

1 

41 

CROD 

152 

1 

42 

CROO 

153 

1 

43 

CROD 

212 

1 

2 

CROD 

221 

1 

12 

CROD 

211 

1 

2 

CROD 

222 

1 

12 

CROO 

231 

1 

22 

CROO 

232 

1 

22 

CROO 

241 

i 

32 

CROD 

242 

1 

32 

CROD 

251 

1 

42 

CROD 

GROSET 

252 

1 

42 

GRID 

1 


.0 

GRID 

2 


.0 

GRID 

3 


.0 

GRID 

11 


40.0 

GRID 

12 


40.0 

GRID 

13 


40,0 

GRID 

21 


80.0 


-30.0 

. 0 

30.0 

-30.0 

.0 

30.0 

-30.0 


.0 

.0 

.0 

.0 

.0 

.0 

.0 


SYNTHESIS 


3456 


30 


GRID 

22 


80. 0 

. 0 

.0 



GRID 

23 


80,0 

30.0 

.0 



GRID 

31 
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. o 



GRID 
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. 0 

.0 



GRID 

33 
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.0 



GRID 

41 
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GRID 

42 
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.0 

.0 



GRID 

43 
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30.0 

.0 



GRID 

51 
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-30.0 

.0 



GRID 

52 


200.0 

. 0 

.0 



GRID 

53 


200.0 

30.0 

.0 



MAT 1 

1 

10. 0*6 


.3 

2.5-4 



PROD 

1 

1 

.3 





EiNlDDAT A 








NASTRAN 

files 

= UMF $ CDC AND 

IBM 




ID 

0EM2032, NASTRAN 






APP D I SP * SUBS 







SOL 2,0 








TIME 3 








CEND 








SUBSTRUCTURE 

PHaSEI 






PASSWORD = MDLSYN 






SOF ( 1 ) 

= FT19 

,500 $ CDC 

AND IBM 




NAME = 

B8ASIC 







SOFPRINT TOC 







ENDSU8S 








title = 

TRUSS 

dynamic ANALYSIS 

USING AUTOMATED MODAL 

SYNTHESIS 


SUBTITLE = NASTRAN DEMONSTRATION PROBLEM 

NO. 2-3-2 



LABEL = 

SUBSTRUCTURE 2, 

RUN 2, 

PHASE 1, 

RF 2 



BEGIN BULK 







CROD 

1 

1 

1 

2 




CROD 

2 

1 

2 
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1 

11 

12 




CROD 

12 

1 

12 
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1 

21 
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1 

22 
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31 

1 

31 
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32 

1 

32 
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41 

1 

41 
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42 

1 

42 
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1 

1 
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CROD 
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1 

2 
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CROD 
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1 

3 

13 




CROD 
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1 

11 

21 




CROD 
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1 

12 

22 




CROD 
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1 

13 

23 




CROD 

131 

1 

21 

31 




CROD 
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1 

22 

32 




CROO 

133 

1 

23 

33 




CROD 
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1 

31 

41 




CROO 

142 

1 

32 

42 




CROD 

143 

1 

33 

43 




CROD 

211 

1 

2 

11 




CROD 

212 

1 

2 

13 




CROD 

221 

1 

12 

21 




CROO 

222 

1 

12 

23 




CROD 

231 

1 

22 

31 




CROD 

232 

1 

22 

33 




CROD 

241 

1 

32 

41 




CROD 

242 

1 

32 

43 




GROSET 
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GRID 
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GRID 
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.0 
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GRID 

3 


-30.0 

.0 
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GRID 
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GRID 

12 
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.0 









31 




GRID 

13 

-30,0 

40.0 

GRID 

21 

30.0 
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GRID 

22 

.0 

80.0 

GRID 

23 

-30.0 

80.0 

GRID 

31 

30.0 

120.0 

GRID 

32 

.0 

120.0 

GRID 

33 

-30.0 

120.0 

GRID 

41 

30.0 

160.0 

GRID 

42 

.0 

160.0 

GRID 

43 

-30.0 

160.0 

MA T 1 

1 

1 0 . 0*6 

.3 

PROD 

ENDDATA 

& 

1 

1 .3 
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.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 


2.5-4 


NASTRAN FILES = UMF $ CDC AND IBM 

ID DEM2033,NASTRAN 

APP DI SP » SUBS 

SOL 3*0 

TIME 5 

CE NO 

SUBSTRUCTURE PHASE2 
PASSWORD = MDL5YN 
SOF(l) = FTI9*500 * CDC AND IBM 
OPTIONS K » M » P 
SOFPRINT TOC 
MREDUCE AdASIC 
NAME MA 
BOUNDARY 5 
FIXED 5 
METHOD 19 
OUTPUT 1»5»6»9»10 
SOFPRINT TOC 
MREDUCE BdASIC 
NAME MB 
BOUNDARY 4 
FIXED 4 
METHOD 29 
OUTPUT 1*5*6*9*10 
SOFPRINT TOC 
MREDUCE CBASIC 
NAME MC 
BOUNDARY 7 
FIXED 7 
METHOD 39 
OUTPUT 1*5,6,9,10 
SOFPRINT TOC 
COMBINE M A » MB , MC 
NAME MCOMB 
TOLERANCE 0.001 
OUTPUT 2*7,12 
COMPONENT MB 
TRANSFORM 20 
COMPONENT MC 
TRANSFORM 40 
SOFPRINT TOC 
MREDUCE MCOMB 
NAME RTRUSS 
BOUNDARY 42 
METHOD 90 
NM AX 18 

OUTPUT 1,5,6*9,10 
SOFPRINT TOC 
ENDSUBS 

TITLE = TRUSS DYNAMIC ANALYSIS USING AUTOMATED MODAL 
SUBTITLE = NASTRAN DEMONSTRATION PROBLEM NO. 2-3-4 
LABEL = MODAL REDUCE* COMBINE* MODAL RECOVERY* RUN 4 
BEGIN BULK 


SYNTHESIS 
PHASE 2* 


RF 3 


32 


ORiGPAL PAQE iS 

OF POOR QUAUTY 


BDYC 

42 

BBASIC 

42 







8DYC 

4 

BBASIC 

3 







BO YC 

5 
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BDYC 

7 

CBASIC 

2 

32 

33 





BOYS 1 

1 

12 

31 





BDYSl 

3 

12 

1 

2 

3 

31 

32 

33 


8DY31 

2 

12 

1 

2 
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BOYS 1 

42 

1 

2 
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EIGR 

1* 

GI V 

.0 


0 



♦ El 

MAX 








♦ E4 

EIGR 

39 

GIV 

.0 

10000.0 


6 



♦ E4 

MAX 








♦ E3 

EIGH 

29 

GIV 

.0 

10000.0 


6 



♦ E3 

MAX 








♦ E2 

EIGR 

90 

GIV 

.0 

10000.0 


20 



♦ E2 

MaX 






30.0 


♦ T2 

TRANS 

20 


240.0 

30.0 

.0 

240.0 

1.0 

+ T2 

200.0 

30.0 

0.0 






♦ T1 

TRANS 

40 


240.0 

.0 

.0 

240.0 

•9 

1.0 

♦ T 1 

ENDDATA 

240.0 

-100.0 

0.0 










NASTRAN 

FILES = UMF $ CDC AND IBM 


ID 

DEM203S. NASTRAN 

■1 

APP DMAP 

* SUBS 


BEGIN 

■ 1 1 1 ■ i ■ 1 1 m i ■ i ii 1 1 1 1 m 1 1 1 1 1 1 1 ■ 1 1 ii 1 1 1 ■ ■ 1 1 iwi 1 1 ■ 


PRECHK 

ALL S 


file 

UDVT =APPEND/TOL= APPEND S 


PARAM 

//*MP Y^/CARDNO/ 0/ 0 $ 


GP1 

GEOMl ♦GE0M2*/GPL*EQEXIN»GPDT*CSTM»8GPDT»SIL/S»N*LUSET/ 5 . N . 


N0GPDT/ALwAY5=-l I 

PLTTRAN 

BGPDT*SIL/BGPDP.SIP/LUSET/S.N.LUSEP s 


PURGE 




COND 

LBL5.N0GPDT S 


GP2 

GE0M2»EQEXIN/ECT $ 


PAR AMI 

PCDB//*PRES*////JUMPPl OT S 


PURGE 

PLTSETX*PLTPAR,GPSETS,ELSETS/JUMPPLOT $ 


COND 

Pi * JUMPPLOT $ 


PLTSET 

pcdb*euexin»ect /pltsetx » pltpar » GPSE TS * ELSETS/S *N*NSIL/S*N* 
jumpplot=-i s 


PRTMSG 

PLTSETX// $ 


PARAM 

//*MPY*/PLTFLG/1/1 s 


PARAM 

//•»MPY*/PFILE/0/0 $ 


COND 

P 1 » JUMPPLOT S 


PLOT 

PLTPAR? GP SETS* EL SETS»CASECC»BGPDT»EQEXIN»SII »»ECT../Pl OTX1/ 


nsil/luset/s»n* jumpplot/s»n*pltflg/s*n»pfile s 

PRTMSG 

PL0TX1// S 


LABEL 

Pi s 


GP3 

GEOM3*EQEXIN»GEOM2/SLT,GPTT/NOGRAV s 


TAl 

ECT »EPT.8GPDT»SIL*GPTT,CSTM/EST*GEI »GPECT» ♦/LUSET/S*N»N0SIMP = 


-l/l/S»N»NOGENt =-l/S*N*GFNEL S 

PURGE 

K4GG»GPST ♦ 0GP5T *MGG*BGG»K4NN*K4FF *K4AA*MNN»MFF *MAA*BNN»BFF » BA A » 


KGGX/NOSIMP/OGPST/GENEL S 

COND 

LBL 1 * NOS I MP S 



33 










ORIGINAL PkGE IS 
OF POOR QUALITY 


p ARAM //»ADO»/NOKGGX/l/Q $ 


PARAM 

PARAM 


PARAM 


tMG 


PURGE 


EMA 


label. 


COMO 


CQND 


GPWG 


OF P 


_EQU_IX 


CO^O 


SMA3 


LABEL 


PARAM 

GP4 


PU RG .EL 


//»ann-»/NQMGG/l/U £ 


//»aod»/nqbgg=-i/i/o S_ 


//»AQO^/NOKAGG/1/0 


'£STfCSTM^MPT»DlT*GE0M2>/KELM>KQICT > MELM » MO I C f » BELM » BO \ C_T/_ 

N *NQKGGX /S»N»NOMGG/S»N»NQBGG/S«M»NOKAG G// C * Y_r CQUPM _ ASS/C » Y _ * _ 

11 w " ~ 1 ' . . . _ — — **...,* ^ * *+ w fr\Tn ? Hi v . 


IN “ 1 1 VJ r\ uu / w t w ^ — — — — : ; . — 

C pg/\p/r t v.rpann/CiYt CPQU AD 1 /C» Y ♦CPQUA02/C»Y »CPTRI A1/C » Yj_ 
CPTRIA^/C«Y»CPTUaE/C*Y < CPQOPLr/CfY »CPTRPLT/C»Y>CPTRBSC_S. 
KGGXtGPST/NOKGGX/MGG/NOMGG i 



GPECT ,KDICT»K£LM/KAGG</NOKAGG 
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JL 

$ 


i flLl «GRDPNT S 


FRR0R3»N0MGG S — — — — — ^fTTTcc 

bgpdp»cstm»eqexin » MGG/OGPWG/ v > Y »grqpn t=-i/c»y» w tm a_s s_ 

OG PWG » ♦ » » »//S»N»CARPNO s 

~LA8E L L6 L1 s — 


_K G GX tKGG/NOG ENL 
fai lltNOGENL S 


£ 


1 P L Ilf tNVVJUTjJ 

GEI >K6GX/KGG/LUSET/N0GENL/N 0SIMP $ 

LBLll $ 


CAS£cLGEOlU»EQEXI N.GPDT,BGPDr»CSTM»GPST/RG»»US£T» ASET/ ^ LU gEIZ. 


S*N. MPCF1/S*N*MPCF2/S*N. SINGLE/S. Nf OMIT/S. N.RtACT/SfN.NSKIP/S, 

n. repeat/s ♦n»noset/s*n»nol/s,n,noa/c,y,asetout/ S,Y,AUTOSPC £ 
OM.GMQ/MPCFl/G Q .Gnn/OMTT/KFSTPST t OP /S T NGL , E - S 



ORK$*AL P&SE IS 
OF POOR QUALITY 





CAS£XX*CSTM>MPT»DIT,EGQYN*SILD* ♦♦BGPOPfTOL* * UPV ♦ EST » XYCDB* 
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APPENDIX A. (concluded) 
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30.0 
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! GRID 
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40.0 
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21 
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i GRTD 

22 
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23 
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I rr f n 

31 


Koni 


WHHMHHHHH 



1 GRID 

32 


.0 

120:0 

.0 



! GRID 

33 


-30.0 


.0 



■m 

41 


30.0 

160.0 

.0 



GRID 

42 


_*0 


.0 



GRID 

43 




_.o 



MAT 1 

1 

10.0+6 


*3 

SBBHHHHHH 


hbbmh 

PROD 

1 

1 

.3 


hb 



DAREA 

980 

42 

1 

1.0*3 




TlO^D x 

/of 

°l^O 




o.o O * ! X. 
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W 
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APPENDIX B. Driver decks and sample bulk data for beam problem of case 1, 


NASTRAN FILES = UMF S CDC AND 

ID DEM2031 » NASTRAN 

APP D I SP » SUBS 

SOL 3,0 

TIME 3 

CEND 

SUBSTRUCTURE PHASE1 

PASSWORD = MDLSYN 

SOF(l) * FT17,600,NEW $ CDC AN[ 

NAME = ABASIC 

SOFPRINT TOC 

ENDSUBS 

TITLE = BEAM DYNAMIC ANALYSIS 


LABEL = 
BEGIN B 
BAROR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
GRDSET 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
PBAR 
MATI 
ENDDATA 
& 


= SUBSTRUCTURE 
BULK 


USING 

PHASE 


10 . 0*6 


0 . 

10 . 

2 °. 

30. 

AO. 

50. 

60. 

70. 

80 . 

90. 

100 . 

.56 


NASTRAN FILES = UMF S CDC ANO IBM 

ID DEM2032, NASTRAN 

APP DISPfSUBS 

SOL 3,0 

TIME 5 

CEND 

SUBSTRUCTURE PHASE2 

PASSWORD = MDLSYN 

SOF(l) = FT17»S00 i CDC AND IBM 

EQUIV ABASIC, B3ASIC 

PRFF T X R 

EQUIV ABASIC, GBASIC 

PREFIX G 

MREDUCE ABASIC 

NAME MA 

BOUNDARY 20 

FIXED 20 

METHOD 1 

MREDUCE B6ASIC 

NAME MB 

BOUNDARY 2 


AUTOMATED 
1, R6 2 

10.0 


MODAL SYNTHESIS 


10.0 


2.591-4 


38 


APPENDIX B. (concluded) 

FIXED 2 
METHOD 2 
MREDUCE GBASIC 
NAME MG 
BOUNDARY 3 
FIXED 3 
METHOD 3 
EQUIV MB » MC 
PREFIX C 
EQUIV MB*MD 
PREFIX D 
EQUIV MB » ME 
PREFIX E 
EQUIV MB *mF 
PREFIX F 

COMBINE MA»M6»MC»MD»M£*MF*MG 

NAME ABCDEFG 

TOLERANCE 0.01 

OUTPUT 2 » 7 » L2 

COMPONENT MB 

TRANSFORM 2 

COMPONENT MC 

TRANSFORM 3 

COMPONENT MD 

TRANSFORM A 

COMPONENT ME 

TRANSFORM b 

COMPONENT MF 

TRANSFORM 6 

COMPONENT MG 

TRANSFORM 7 

MREDUCE ABCDEFG 

NAME BEAM 

BOUNDARY 20 

METHOD 22 

OUTPUT 1*5*6,9,10 

SOFPRINT TOC 

ENDSUBS 

TITl E=BEAM DYNAMIC ANALYSIS USING AUTOMATED MODAL SYNTHESIS 
SUBTiTLE=NASTRAN DEMONSTRATION PROBLEM NO. 2-3-2 

label=modal — " * — * 

BEGIN BU^K 


* , nni* i i wti r nyuLW." nut 

REDUCE* COMB I NE, MODAL, RECOVER Y » RUN2 * PHASE2 


BDYC 
BDYC 
BDYC 
BDYS1 
BDYS1 
BDYS1 
TRANS 

♦ T2 
TRANS 

♦ T3 
TRANS 

♦ TA 
TRANS 

♦ T5 
TRANS 

♦ T6 
TRANS 

♦ T7 
EIGR 

♦ El 
EIGR 

♦ E2 
EIGR 

♦ E3- 
EIGR 

♦ E22 
ENDOATA 

♦ 


2 

3 

20 

30 

A0 

50 

2 

150. 

3 

250. 

A 

350. 

A50 . 
6 

550. 

650 . 
1 

MAX 

2 

MAX 

3 

MAX 

22 

MAX 


BBASIC 

GBASIC 

ABASIC 

126 

126 

126 

0. 
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APPENDIX C. Driver decks and sample bulk data for beam problem of case 2. 


NASTRAN FILES = UMF % CDC AND IBM 

ID DEM2031 .NASTRAN 

APP D I SP » SUBS 

SOL 3*0 

TIME 3 

CEND 

SUBSTRUCTURE PHASE1 
PASSWORD = MDLSYN 
S0F(1) = FT 1 7.500. NEW $ CDC 
NAME = ABASIC 
SOFPRINT TOC 
ENDSU8S 

title = beam dynamic analysis 
■“ = SUBSTRUCTURE 1* RUN I* 

BULK 


LABEL 
BEGIN 
BAROR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
GROSET 
GRID 
GRID 
GRID 
GRID 
GRID 
OHIO 
GRID 
GRID 
GRID 
GRID 
GRID 
PBAR 
MAT1 
ENDDATA 
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NASTRAN FILES = UMF S CDC AN0 IBM 

ID DE M20 32 » NASTRAN 

APP D I SP * SUBS 

SOL 3.0 

TIME 5 

CEND 

SUBSTRUCTURE PHASE2 

PASSWORD = MDLSYN 

SOF(l) = FT17.500 S CDC AND IBM 

EQU I V ABASIC.9BASIC 

PREFIX B 

EQUIV ABASIC.CBASIC 
PREFIX C 

EQUIV ABASIC .DBAS I C 
PREFIX D 

EQUIV ABASIC. EBASIC 
PREFIX E 

EQUIV ABASIC. F3ASIC 
PREFIX F 

EQUIV ABASIC.GBASIC 
PREFIX G 
MREDUCE ABASIC 
NAME MA 


AUTOMATED 
1. Rb 2 

10.0 


MODAL SYNTHESIS 


10.0 


0.0 


345 


0. 
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BOUNOARY 20 
FIXED 2 0 
METHOD 1 
MREDUCE BBASIC 
iM A M E M ti 
BOUNDARY 2 
FI XED 2 
METHOD 2 
MREDUCE CBASIC 
NAME MC 
BOUNDARY 7 

FIXED 7 
METHOD 2 
MREDUCE DBAS I C 
NAME MD 
BOUNDARY b 
FIXED 8 
METHOD 2 
MREDUCE EBASIC 
NAME ME 
BOUNDARY 9 
FIXED 9 
METHOD 2 
MREDUCE FBASIC 
NAME MF 
BOUNDARY 11 
FIXED 11 
METHOD 2 
MREDUCE GBASIC 
NAME MG 
BOUNDARY 3 
FIXED 3 
METHOD 3 
COMBINE MA »M8 
NAME AB 
TOLERANCE 0.01 
OUTPUT 2,7,12 
COMPONENT MB 
TRANSFORM 2 
MREDUCE AB 
NAME MAB 
BOUNDARY 10 
FIXED 10 
METHOD 22 
COMBINE MAB , MC 
NAME ABC 
TOLERANCE 0.01 
OUTPUT 2,7,12 
COMPONENT MC 
TRANSFORM 3 
MREDUCE ABC 
NAME MABC 
BOUNDARY 21 
FIXED 21 
METHOD 22 
COMBINE MABC *MD 
NAME ABCD 
TOLERANCE 0.01 
OUTPUT 2,7,12 
COMPONENT MD 
TRANSFORM A 
MREDUCE ABCD 
NAME M ABCD 
BOUNDARY 22 
FIXED 22 
METHOD 22 


°ki&nal paez is 

0F POOR quality 
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combine mabcd»me 
name A8CDE 
TOLERANCE 0.Q1 
OUTPUT 2*7*12 
COMPONENT ME 
TRANSFORM 5 
MREDUCE ABODE 
NAME MABCDE 
BOUNDARY 23 
FIXED 23 
METHOD 22 
COMBINE MABCDE*MF 
NAME ABCDEF 
TOLERANCE 0.01 
OUTPUT 2*7,12 
COMPONENT MF 
TRANSFORM 6 
MREDUCE ABCDEF 
NAME MABCDEF 
BOUNDARY 24 
FIXED 24 
METHOD 22 

COMBINE MaBCDEF *MG 
NAME ABCDEFG 
TOLERANCE 0.01 
OUTPUT 2*7*12 
COMPONENT MG 
TRANSFORM 7 
MREDUCE ABCDEFG 
NAME BEAM 
BOUNDARY 20 
METHOD 22 
OUTPUT 1*5*6,9,10 
SOFPRINT TOC 
ENOSUBS 

I, I l I^? 8 l AM .9Xf!!^ MI ^ ANA l-YSIS USING AUTOMATED mooal synthesis 
SU8T I TlE=NASTR AN DEMONSTRATION PROBLEM NO. 2-3-2 
LABEL=MODAL REDUCE*C0MBINE*M0DAL»REC0V5RY»RUN2»Pa S52 
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APPENDIX C. (concluded) 
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APPENDIX D. Driver decks and sample bulk data fro beam problem of case 3. 


NASTRAN FILES = UMF S CDC AND IBM 

ID DEM2031 f NASTRAN 

APP Q I SP * SUBS 

SOL 3 » 0 

TIME 3 

CEND 

SUBSTRUCTURE PHASE1 
PASSWORD = MDLSYN 

SOF { 1 ) = FT 17 f 50 Of NEW S CDC AND IBM 
NAME = ABASIC 
SOFPRINT TOC 
ENDSUBS 

TITLE = BEAM DYNAMIC ANALYSIS 


LABEL 
BEGIN 
BAROR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
CBAR 
GRDSET 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
GRID 
PBAR 
MAT1 
ENDDATA 
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automated 

If R6 2 

10.0 


MODAL 


10.0 


NASTRAN FILES = UMF S CDC AND IBM 

ID DEM2032f NASTRAN 

APP OISP f SUBS 

SOL 3 f 0 

TIME 10 

CEND 

SUBSTRUCTURE PHASE2 

PASSWORD = MDLSYN 

SOF ( 1 ) = FT 17 *500 $ CDC AND IBM 

EQUIV ABASICfBBASIC 

PREF IX B 

EQUIV ABASICfCBASIC 
PREFIX C 

EQUIV ABASICfDBASIC 
PREFIX D 

EQUIV ABASICfEBASIC 
PREFIX E 

EQUIV ABASICf FBASIC 
PREFIX F 

EQUIV ABASICfGBASIC 
PREFIX G 
MREDUCE ABASIC 


0. 

0. 

0. 

0. 

0. 

0. 

0 . 

0 . 

0 . 

0 . 
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SYNTHESIS 

0.0 1 


34d 


44 


NAME MA 
BOUNDARY BO 
FIXED 20 
METHOD 1 
MREDUCE BBASIC 
NAME MB 
BOUNDARY 2 
FIXED 2 
METHOD 2 
MREDUCE CbASIC 
NAME MC 
BOUNDARY 7 
FIXED 7 
METHOD 2 
MREDUCE DBASIC 
NAME MD 
BOUNDARY B 
FIXED « 

METHOD 2 
MREDUCE EBASIC 
NAME ME 
0OUNOARY 9 
FIXED 9 
METHOD 2 
MREDUCE FBASIC 
NAME MF 
BOUNDARY 11 
FIXED 11 
METHOD 2 _ 

MREDUCE GBASIC 
NAME MG 
BOUNDARY 3 
FIXED 3 
METHOD 3 
COMBINE MA »MB 
NAME AB 

TOLERANCE, 0.01 
OUTPUT 2 » 7 * 12 
COMPONENT MB 
TRANSFORM 2 
MREDUCE AB 
NAME MAS 
BOUNDARY 10 
FIXED 10 
METHOO 22 
COMBINE MCtMD 
NAME CD 

tolerance 0.01 
OUTPUT 2 * 7 » 12 
COMPONENT MD 
TRANSFORM 2 
MREDUCE CD 
NAME MCD 
BOUNDARY 15 
FIXED 15 
METHOD 22 
COMBINE ME *MF 
NAME EF 

TOLERANCE 0.01 
OUTPUT 2 » 7 * 12 
COMPONENT MF 
TRANSFORM 2 
MREDUCE EF 
NAME MEF 
BOUNDARY 25 
FIXED 25 
METHOO 22 
COMBINE MAB »MCO 


NAME ABCD 
TOLERANCE 0,01 
OUTPUT 2,7*12 
COMPONENT MCO 
TRANSFORM 3 
MREDUCE AtJCD 
NAME MA0CD 
BOUNDARY 30 
FIXED 30 
METHOD 25 
COMBINE MEF *M6 
NAME EF6 

tolerance 0.01 

OUTPUT 2*7,12 
COMPONENT MG 
TRANSFORM 3 
MREDUCE EFG 
NAME MEFG 
BOUNDARY 35 
FIXED 35 
method 25 

COMBINE mabcd*mefg 
name abcdefg 
tolerance 0.01 
OUTPUT 2,7,12 
componet mefg 
TRANSFORM 5 
MREDUCE ABCDEFG 
NAME BEAM 
BOUNDARY 20 
METHOD 25 
OUTPUT 1,5,6,9,10 
SOFPRINT TOC 
ENDSUBS 

TIT!E=BEAM OYNAMIC ANALYSIS USING AUTOMATED MODAL SYNTHESIS 
SUBTITlE=NASTRAN DEMONSTRATION PROBLEM NO. 2-3-2 
LABEL=MODAL REDUCE , COMBINE, MODAL , RECOVERY, RUN2,PHASE2 
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APPENDIX D. (concluded) 
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